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ABSTRACT 

A refined model for the unusual Type Ic supernova 1998bw, discovered as the optical counterpart of 
GRB980425, is presented, and synthetic Ught curves and spectra are compared with the observations. 
The first 30 days of the light curve and the broad line features of the spectra can be reproduced with the 
hydrodynamical model of the explosion of a 14M0 C+0 star, the core of a star with initial mass 4OM0, 
assuming that the explosion was very energetic (kinetic energy Ek = 5 x 10^^ erg) and that QAMq of 
^^Ni were synthesized. At late times, however, the observed light curve tail declines more slowly than 
this energetic model, and is in better agreement with a less energetic (Sk = 7 x 10^^ erg) one. This 
shift to a less energetic model may imply that the inner part of the ejecta has higher density and lower 
velocities than the model with Ey^ = 5 x 10^^ erg, so that 7-rays deposit more efficiently. An aspherical 
explosion can produce such a structure of the ejecta. We also study detailed nucleosynthesis calculations 
for hyper- energetic supernova explosions and compare the yields with those of normal supernovae. 

Subject headings: gamma-ray burst — supernovae — supernovae: individual (SN1998bw) — 
nucleosynthesis 



1. INTRODUCTION 

On 1998 April 25, GRB980425 triggered the Narrow 
Field Instrument (NFI) and Wide Field Camera (WFC) 
detectors on board BeppoSAX (Soffita et al. 1998). 
SN 1998bw was discovered within the WFC error box in 
the optical (Calama et al. 1998a) and radio wavelength 
bands (Kulkarni et al. 1998a) only 0.9 and 3 days after 
the date of the GRB, respectively. The X-ray afterglow 
detected by the NFI within the error box of CRB980425 
is consistent with SN 1998bw (Pian et al. 1999). The 
small likelihood of finding a supernova and a CRB in such 
a small field over such a short interval of time suggests 
that SN 1998bw and GRB980425 are related (Galama et 
al. 1998b; Kulkarni et al. 1998b). 

The early optical spectra of SN 1998bw apparently lack 
dominant line features, displaying only some broad emis- 
sions at 500, 620 and 800 nm (Galama et al. 1998b; 
Stathakis et al. 2000; Patat et al. 2000). These were 
shown to be the emission components of P-Cygni profiles 
of very broad line blends. They are caused mostly by 
Fell lines in the blue, by Sill near 600 nm and by 01 -I- 
Call near 720 nm (Iwamoto et al. 1998, hereafter IMN98). 
The absence of any hydrogen line after subtraction of the 
galaxy background and the fact that the Si II line at 615 
nm has a very large velocity indicate that the supernova 
is neither a Type II nor a normal Type la (Galama et al. 
1998b). Following the conventional classification scheme 
(Filippenko 1997), the lack of strong He I features led 
Patat & Piemonte (1998) to conclude that SN 1998bw is 
likely to be a Type Ic supernova (SN Ic) rather than a 
Type lb. 

The fight curve of SN 1998bw (Galama et al. 1998b) 



showed a very early rise and reached a peak at '^ 17 days 
(in the V band) after the explosion, and then has been 
declining exponentially with time (McKenzie & Schaefer 
1999; Patat et al. 2000). This clearly indicates that 
the light curve is not a typical optical afterglow of a 
gamma-ray burst, but it is powered by the radioactive de- 
cay of ^^Ni and ^^Co as in usual supernovae (IMN98). 
The distance modulus to SN 1998bw is estimated as /x = 
32.89 mag, so that the peak luminosity of SN 1998bw is 
^ 1 X 10^^ ergs sec~^, which is about ten times brighter 
than typical SNe Ib/Ic (Clocchiatti & Wheeler 1997). 
To achieve such a high luminosity, a large amount of 
^^Ni must have been synthesized in SN1998bw (IMN98; 
Woosley, Eastman & Schmidt 1999), again about ten times 
as much as that produced in typical core-collapse-induced 
supernovae. The very broad spectral features and the light 
curve shape led various groups to the conclusion that SN 
1998bw had a very large kinetic energy of explosion Ek 
(IMN98; Woosley et al. 1999; Branch 2000). 

IMN98 constructed models of the core-collapse-induced 
explosion of C-l-0 cores of initially massive stars that had 
lost their hydrogen and helium-rich layers before the ex- 
plosion. Among those models, the energetic explosion 
{Ek ^ (2 - 5) X 10^2 gj.gs) ^f ^ C-hO star of 13.8 Mq 
successfully fit the first 60 days of the light curve and 
the early spectra, although the synthetic spectral lines 
were still narrower than the very broad observed features 
(IMN98). Since the kinetic energy is more than one order 
of magnitude larger than the energy of typical supernovae, 
SN 1998bw was caUed a "hypernova" (IMN98), a term we 
use to describe events with Ek > 10^^ erg without specify- 
ing whether the central engine is a coUapser (MacFadyen 
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& Woosley 1999), a magnetar (Nakamura 1998; Wheeler 
et al. 2000) or a pair-instability supernova. 

Interestingly, photometry after ^ 60 days showed that 
SN 1998bw declined significantly more slowly than the rate 
predicted by the model of IMN98 (McKenzie & Schaefer 
1999; Patat et al. 2000). Also the bolometric light curve 
has been constructed up to day 500 (Patat et al. 2000; 
SoUerman et al. 2000). Therefore, we have recomputed 
the light curve of SN 1998bw using progenitors of different 
masses and explosions of different energies, and provide 
better estimates for these parameters. 

The hydrodynamical models for SN 1998bw are de- 
scribed in §0. Sections ^ and H are devoted to light curve 
models and synthetic spectra models, respectively. Ex- 
plosive nucleosynthesis in SN 1998bw is discussed in §|[ 
Finally, the nature of this peculiar supernova is summa- 
rized in §0, with the emphasis being placed on possible 
evidence that the explosion was aspherical. Preliminary 
results have already been reported by Nakamura et al. 
(1999a, 2000) and Nomoto et al. (2000ab). 

2. HYDRODYNAMICAL MODELS FOR SN 1998BW 

Hydrodynamical models are constructed as follows. 
C-l-0 stars are chosen as progenitors as in IMN98. Light 
curves and spectra are computed for various C-l-0 star 
models with different values of the kinetic energy Ek and 
the ejecta mass Mgj. These parameters can be constrained 
by comparing calculated light curves, synthetic spectra, 
and photospheric velocities with the observational data of 
SN 1998bw. 

(1) In the ordinary, low energy SN Ic model (model 
C0138E1), a C-l-0 star with a mass Mco = IS.SM© 
(which is the core of a 40 AIq main-sequence star; 
Nomoto & Hashimoto 1988; Nomoto et al. 1997) 
explodes with Ek — 1.0 x 10^^ ergs and Mej = 
Mco - Mcut ^ 12Mq. Mcut (= 2 Mq in this case) 
denotes the position of the mass cut, which corre- 
sponds to the mass of the compact star remnant. 
This is either a neutron star or a black hole, depend- 
ing on Mcut- 

(2) For the hypernova models CO138E50, CO138E30, 
and C0138E7, the progenitor C-f O star of Mco = 
I3.8M0 is the same as in C0138E1. These models 
have different explosion energies: Ek — 5 x 10^^ 
erg (CO138E50), 3 x 10^^ gj.g (CO138E30), and 
7 X 10^1 erg (C0138E7). The ejecta mass is Mej ~ 
10 - II.5M0, i.e., Mcut ~ 2.5 - 4 Mq. The param- 
eters of the models are summarized in Table H. The 
position of the mass cut is chosen so that the ejected 
mass of ^^Ni is the value required to explain the ob- 
served peak brightness of SN 1998bw by radioactive 
decay heating. The compact remnant in these hyper- 
nova models may well be a black hole, because Mcut 
can exceed, sometimes significantly, the maximum 
mass of a stable neutron star. 

The hydrodynamics at early phases was calculated using 
a Lagrangian PPM code (Colella & Woodward 1984). All 
models are spherically symmetric. The explosion is trig- 
gered by depositing thermal energy in a couple of zones 
just below the mass cut so that the final kinetic energy 



has the required value. A strong shock wave forms and 
propagates toward the surface. Explosive nucleosynthesis 
takes place behind the shock wave. Radioactive ^^Ni is 
produced in the deep, low velocity layers of the ejecta. 

In SN II and SN lb models, it has been demonstrated 
that Rayleigh- Taylor instabilities develop at the H/He and 
He/ C-l-0 interfaces and induce mixing of elements in the 
ejecta (e.g. Arnett et al. 1989; Hachisu et al. 1991, 1994; 
Iwamoto et al. 1993). Although bare C-l-0 stars, the 
progenitors of SNe Ic, lack composition interfaces with 
such pronounced density jumps, Rayleigh- Taylor instabil- 
ities can be driven by neutrino heating and develop at the 
Ni-l-Si/0 interface (Kifonidis et al. 2000). In addition, 
polarization measurements suggested that SNe Ic, and hy- 
pernovae in particular, are asymmetric explosions (Wang 
et al. 1999; Danziger et al. 1999; Patat et al. 2000). The 
asymmetry of the explosion may reinforce the instability 
and bring heavy elements up to high velocity layers. This 
is particularly true of jet-like explosions, where it is likely 
that extensive mixing takes place in velocity space. In view 
of the large uncertainties in our knowledge of the mixing 
process, we assume that the material ejected is uniformly 
mixed out to a velocity v = Umix- 

The hydrodynamical models become homologous {v oc 
r) at i ~ 250 sec, and are then used as input for a radia- 
tion transfer code. The lines in Figure |l| show the density- 
velocity distribution of the homologously expanding ejecta 
(top) and the enclosed mass Air as a function of velocity 
(bottom) for models CO138E50 (solid), CO138E30 (long- 
dashed), C0138E7 (short-dashed), and C0138E1 (dash- 
dotted). 

3. LIGHT CURVE MODELS 

Synthetic light curves are computed with a radiative 
transfer code (Iwamoto et al. 2000) which takes into ac- 
count the balance between photo-ionizations and recombi- 
nations and includes a simplified treatment of line opacity. 
The width of the light curve peak, tlc, depends on i?K and 
Mej approximately as tlc ^ {k/cY/'^M^^ E^ ' , where k 
and c are the optical opacity and the speed of light, re- 
spectively (Arnett 1982). This means that the light curve 
can be reproduced with different explosion models that 
have the same values of M^-E^^. However, these param- 
eters can be further constrained from both photospheric 
velocities and spectra because the velocity scales roughly 
as Al- ' E-jl . The light curve shape depends also on 
the distribution of the radioactive heating source ^^Ni, for 
which we examine the dependence on Umix (§0)- 

Our calculations are compared with the bolometric light 
curve of SN 1998bw constructed by Patat et al. (2000), 
which uses a redshift distance n = 32.89 mag {d ~ 37.8 
Mpc, Ho = 65 km s^^ Mpc^^) and Ay = 0.05 mag. The 
time of core collapse is set at the detection of GRB 980425. 
Figure ^ shows the bolometric light curves of the energetic 
models CO138E50 (solid) and CO138E30 (long-dashed), 
and Figure || the less energetic models C0138E7 (short- 
dashed) and C0138E1 (dash-dotted). 

Figure ^ shows the evolution of the calculated photo- 
spheric velocities of CO138E50 (sohd line), CO138E30 
(long-dashed line), C0138E7 (short-dashed line), and 
C0138E1 (dash-dotted line) compared with those ob- 
tained from spectral models (filled circles), and with the 
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observed velocity of the Si II 635.5 nm doublet measured 
in the spectra at the absorption core (open circles; Patat et 
al. 2000), and that of the Ca II H + K doublet measured i n 
the spectrum of May 23 (square; Patat fc Picmontc 1998| ). 



3.1. Early Phase 

The early part of the light curve (t < 25 days) is well 
reproduced by the two energetic models, CO138E50 and 
CO138E30. The difference between the two curves is too 
small to see on the scale of Figure because the light curve 
width scales as Ej^ . The early part of the light curve 
also depends on the ^^Ni distribution in the ejecta. To 
reproduce the early sharp rise of the light curve, extensive 
mixing of ^^^^^ ^^^ ^^ „^^.^ ^ 22, 000 km s^^ for CO138E50 
and Wmix = 30, 000 km s^^ for CO138E30, is required. A 
better fit can be obtained if we adopt the following ad hoc 
S'^Ni distribution: X(56Ni) = 0.79 at •;; < 11,000 km s-\ 
0.011 at 1) = 11 - 17,000 km s~\ and 0.032 at u = 17 - 
40,000 km s~^, where X(^^Ni) denotes the mass fraction of 
^^Ni. The result is shown in Figure R. Such a non-uniform 
distribution of ^^Ni might be due to the mixing caused by 
Rayleigh- Taylor instability (§^), or reflect a complicated 
structure of a possible jet-like ejecta (§^. The ^^Ni mass 
is determined to be 0.4 Mq from the fitting of the max- 
imum luminosity and the date of maximum. (Note that 
IMN98 slightly overestimated the mass of ^^Ni because 
they assumed that the bolometric correction is negligible 
and adopted a different absorption and distance.) 

In contrast, the early light curves of the less energetic 
models (C0138E7 and C0138E1) evolve too slowly, reach- 
ing the maximum too late as compared with the obser- 
vations (Figure 0), even if ^^Ni is distributed uniformly 
throughout the ejecta. 

The photospheric velocities Vph provide clearer diag- 
nostics to distinguish between CO138E50 and CO138E30. 
In particular, in the early phase (i < 20 days), Wph of 
CO138E50 is in good agreement with the observed Wph, 
while Vph of CO138E30 is clearly too low to be consistent 
with the observations (Fig. |4|). 

In the early phase, therefore, model CO138E50 shows 
the best agreement with both the light curve and the pho- 
tospheric velocities of SN 1998bw. The synthetic spectra 
also require the most energetic model CO138E50 as will 
be shown in pl CO138E50 has the same mass as the best 
model in IMN98 but a larger Ek, which is necessary espe- 
cially to improve the fit to the spectra (§Q). 

3.2. Intermediate Phase 

At intermediate phases (t ~ 25 - 200 days), the decline 
rate of the light curve is determined mainly by the fraction 
of the 7-rays emitted by ^^Co decay which are trapped 
in the ejecta. The optical depth of the ejecta to the 7- 
rays scales as n^pR ex MR~^ oc M'^E^^t~^. Thus the 
behaviour of the various models can be seen more easily 
than at earlier phases as follows: 

(1) The light curve of CO138E50 is consistent with 
SN1998bw until day 50 but declines faster than 
the observation afterwards (Figure 0; McKenzie & 
Schaefer 1999; Patat et al. 2000). The light 
curve of CO138E30 declines more slowly than that 
of CO138E50, and is in better agreement with 



SN1998bw, but it still dechnes faster than the ob- 
servations. 

(2) The photospheric velocity Uph shows a similar ten- 
dency. Figure § shows that Vph of CO138E50 is the 
best fit to the data for the first 20 days, but after- 
wards Wph of CO138E30 gives as good a fit. The 
other two models do not fit Vph at all. 

(3) Figure shows that the apparently exponential de- 
chne of SN1998bw after day 60 (McKenzie & Schae- 
fer 1999; Patat et al. 2000) is well reproduced by 
the lower energy model C0138E7. In this model, 
7-ray trapping is more efficient than CO138E50 and 
CO138E30, but the observed low flux level at this 
phase requires a reduced ^^Ni mass of 0.28Mo (dot- 
ted hne). Note, however, that 0.28Mo of ^"^Ni is too 
small to reproduce the observed light curve maxi- 
mum. Note also that because of the interplay of Ek 
and M(Ni), it is difficult to establish those parame- 
ters from the intermediate phase light curve alone. 

(4) The normal SN Ic model C0138E1 in Figure | is 
slow enough to trap most of the 7-rays emitted from 
the ^^Co decay. The decline of the light curve is 
therefore too slow compared with the observed rate. 

These comparisons between SN 1998bw and the model 
light curves and Wph in Figures 0, 0, & indicate that 7- 
ray deposition after about day 50 in SN 1998bw is more 
efficient than predicted by CO138E50 and E30, although 
these models give an appropriate description of the early 
light curve and spectra. Higher deposition can be achieved 
if there exists a significant amount of low-velocity and 
high-density material. Indications of the presence of a 
low-velocity, high-density region suggest that the ejecta 
distribution is not spherically symmetric as will be dis- 
cussed in §0. We note that such indications are also found 
for SN 1997ef (Iwamoto et al. 2000; Mazzah, Iwamoto, & 
Nomoto 2000a), a lower-energy analogue of SN 1998bw. 

3.3. Late Phase 

After day ~ 200 the decline of the model light curve be- 
comes slower, and it approaches the half-life of ^^Co decay 
around day 400. At t > 400 days most 7-rays escape from 
the ejecta. The 7-ray deposition fraction at 400 days is 1%, 
1.5%, 6%, and 13% for CO138E50, E30, E7, and El, re- 
spectively. On the other hand, kinetic energies of positrons 
emitted from the decay of ^^Co are supposed to be fully 
thermalized because of the postulated weak magnetic field 
(e.g., Colgate & Petchek 1979). Therefore, positron kinetic 
energy deposition determines the light curve at t > 400 
days (dotted line in Figure 0). Here the luminosity by 
positron deposition is given as 



L(Co,e+) = 1.4 X lO^'^erg/s (A/(Ni)/M0) 
xexp(-t/111.26(day)) x 0.035, 



(1) 



where the positron fraction in energy is 3.5 percent (e.g., 
Axelrod 1980). (See Cappellaro et al. 1997 and Milne et 
al. 1999 for the light curves including positron escape.) If 
the observed tail should follow the positron-powered light 
curve, the ^^Co mass could be determined directly. Since 
positron deposition should occur almost on the spot, this 
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determination does not depend much on any asphericity 
of the ejecta. 

On June 11, 2000, which corresponds to a SN epoch 
oi t = 778 days, HST observations detected a point-like 
source at the position of SN1998bw (Fynbo et al. 2000). 
The observed magnitude (V-mag 25.41 ± 0.25) is con- 
sistent with the prediction of C0138E7 (Figure ||) but 
brighter than CO138E50 (Figure ||). The observed lumi- 
nosity would even be higher than C0138E7, if the bolo- 
metric correction was significant. In any case, these com- 
parisons suggest that the positron contribution is not dom- 
inant yet around day 800. 

4. SYNTHETIC SPECTRA 

In Figure we show the synthetic spectra obtained for 
the same 3 epochs fitted in IMN98 (continuous lines) com- 
pared with spectra observed at ESO (bold lines). We used 
model CO138E50 and computed synthetic spectra with a 
Monte Carlo model (Mazzali & Lucy 1993), improved with 
the inclusion of photon branching and a new extended and 
improved line list (Lucy 1999; MazzaH 2000). The syn- 
thetic spectra were computed using the luminosity derived 
from the light curve, a distance modulus of /x = 32.89 mag 
and Ay = 0.05. The assumption of low reddening is sup- 
ported by the upper limit of 0.1 A in the equivalent width 
of the Na i D line obtained from high-resolution spectra 
(Patat et al. 2000). The observed spectra used here (solid 
Hues) are the 'definitive', fully reduced version of the same 
ESO spectra shown in IMN98, and are calibrated with re- 
spect to the V photometry. The residual correction factors 
for the other bands are usually very close to 1 , but they are 
^ 1.1 for B in the May 11 and 23 spectra. Therefore, the 
new models (i.e. 3 epochs of CO138E50) have somewhat 
different parameters than those of IMN98. Both the lumi- 
nosities and the photosphcric velocities are larger than in 
IMN98. The photosphcric velocity is now in better agree- 
ment with the measured velocity of the Si ii line (IMN98, 
Figure 3). 

The synthetic spectra clearly improve over those of 
IMN98, especially at the earliest epochs. Absorptions not 
caused by broad blends of many lines of moderate strength, 
such as the Si ii feature near 6OOOA and, in particular, the 
O I-|-Ca II feature between 7000 and 8OOOA, are now much 
broader, in significantly better agreement with the data. 
Nevertheless, the blue sides of those absorptions are still 
too narrow, indicating that even the new model CO138E50 
may not contain enough mass at the highest velocities. 

Therefore we introduced an arbitrary change to the 
original CO138E50 density structure. Several possibili- 
ties were tested, and improved results were obtained when 
the density slope was reduced from p oc r^^ to p oc r^^ 
at w > 30,000 km s~^. This does not introduce a sig- 
nificant change in M^j, and increases Ek by only about 
10%, but it does increase the density at high velocities, 
leading to significant absorption at w ~ 60,000 km s~^ in 
the strongest lines, especially the Ca 11 IR triplet, extend- 
ing the absorption troughs to the blue. The corresponding 
synthetic spectra are shown as the dotted lines in Figure pi 
The effect of the change is of course largest at the earliest 
epochs. Although the overall agreement with the observed 
spectra is better, several problems remain, the most severe 
of which is clearly the excessive strength of the O i line at 
7200 A on May 11 and 23. The composition of the highest 



velocity ejecta is dominated by O in our ID models, and it 
is difficult to make that line become weaker. On 23 May, 
the synthetic Ca il IR triplet matches the weak feature at 
8OOOA, which is first seen on 11 May and which contin- 
ues to grow until it finally causes the wavelength of the 
absorption minimum of the entire broad feature to shift 
to - 8200A (Patat et al. 2000). This is rather a peculiar 
behavior, because on 3 May the O i and Ca 11 lines had to 
blend much more to give rise to the observed broad fea- 
ture, which then had a minimum at 7000A. The synthetic 
O I line is too strong and too fast. The core of the broad 
feature, if it is interpreted as O i 7774A, indicates a ve- 
locity of 10500 km s"^ on May 11 and 6000 km s^^ on 
May 23. This is significantly lower than the corresponding 
velocity of the model photosphere. Actually, the defini- 
tion of a photosphere at very red wavelengths is not very 
accurate, because the density of spectral lines is low and 
so line opacity does not define a pseudo-continuum in that 
region. Therefore, the observations probably indicate that 
a large fraction of the O is located at low velocities. 

A very fiat (p oc r~^) density distribution was also used 
by Branch (2000) to fit the spectrum of SN 1998bw. This 
dependence is however too flat when we use our MC model, 
because the ionization of, e.g., Ca 11 does not fall as steeply 
as he assumed. On the other hand. Branch's value of Ek 
(5 X 10^^ erg) is similar to ours, but he quotes a mass of 
6M0 above 7000 km s^^, while in our case the mass above 
that velocity is as large as ~ IOMq. Such a fiat density 
distribution at high velocities is also required to fit the 
spectrum of another hypernova, SN 1997ef (Mazzali et al. 
2000a). This might indicate that the progenitors of these 
hypernovae underwent very extensive mass loss; the outer 
density structure of the ejecta is fiatter than that of ordi- 
nary giants, and comparable to that of a mass-losing star. 

Clearly, a definitive solution has not been found yet. 
It is quite possible that only by taking into account de- 
partures from spherical symmetry will it be possible to 
obtain a really accurate fit to the spectra. Nevertheless, 
considering the complexity of the problem, our fits at least 
demonstrate that a large Ek is necessary and that the O- 
dominated composition of the SN Ic model yields quite a 
reasonable reproduction of the observations. 

5. EXPLOSIVE NUCLEOSYNTHESIS 

We calculated explosive nucleosynthesis using a de- 
tailed nuclear reaction network (Thielemann, Nomoto & 
Hashimoto 1996; Nakamura et al. 1999b). Our calcu- 
lations are performed in two steps. The first step is a 
hydrodynamical simulation of the explosion with a small 
nuclear reaction network containing only 13 alpha nuclei 
(^He, 12c, 160, 20Ne, 24Mg, 2831, ^2^^ 36^^, 40ca, 44^1^ 
4^Cr, ^2pe, and ^^Nj), as described in §||. In the second 
step, post-processing calculations are performed at each 
mesh point of the hydrodynamical model with the ex- 
tended reaction network (Hix & Thielemann 1996), which 
contains 211 isotopes up to ^^Ge. 

The top panel of Figure R shows the isotopic composi- 
tion of the ejecta of the hypernova model CO138E50 as 
a function of the enclosed mass Mr and of the expansion 
velocity. The nucleosynthesis in other hypernova and su- 
pernova models (CO138E30, C0138E7, C0138E1) is also 
shown in Figure M for comparison. The yields of the hy- 
pernova and supernova models are summarized in Table 0. 
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Table || and Figure || give respectively more detailed yields 
and the abundances of stable isotopes relative to the so- 
lar values for model CO138E50. From these figures and 
tables, we note the following characteristics of nucleosyn- 
thesis of hyper-energetic explosions compared with normal 
energy explosions. 

(1) The complete Si-burning region where ^^Ni is pro- 
duced is extended further out (in mass coordinates) 
as the explosion energy increases. How much pro- 
cessed matter is ejected from this region depends 
on the mass cut. Compared with normal core- 
collapse supernovae, the much larger amount of ^^Ni 
(~ OAMq) observed in SN1998bw impUes that the 
mass cut is deeper, so that the elements synthesized 
in this region, such as ^^Cu, ^^Zn, and ^^Ge (which 
decay into ^^Co, ^■^Cu, and ^^Zn, respectively), are 
ejected more abundantly. Among hypernova models, 
more energetic models produce more ^^Ni in the in- 
complete Si burning region (See (3) below), and thus 
Mcut is larger to eject ^ OAMq ^^Ni as constrained 
from the light curve of SN 1998bw (Figure 0). 

(2) In the complete Si-burning region of the hypernova 
models (CO138E50/E30), elements produced by a- 
rich freezeout are enhanced because nucleosynthesis 
proceeds at lower densities than in C0138E1. Fig- 
ure ^ clearly shows a trend that a larger amount 
of ''He is left in more energetic explosions. Hence, 
the mass fractions of the species synthesized through 
a-particle capture, such as '*''Ti and ^*Cr (which 
decay into '*''Ca and "^^Ti, respectively) is larger 
in CO138E50/E30/E7 than C0138E1. The inte- 
grated mass of these species depends on Mcut- For 
CO138E50, the ejected mass of *'*Ti is smaller than 
in other models because of its larger Mcut- Note that 
the ^He produced even in the most energetic mod- 
els has a velocity significantly smaller than that of 
the He shell identified in SN 1998bw at a velocity of 
18300km s^i by Patat et al. (2000) on the basis of 
near-IR spectra. Spectral evidence for low velocity 
He in SN 1998bw is unclear, however the distribution 
of He depends on mixing in velocity space. 

(3) The incomplete Si-burning region is more massive 
in more energetic explosions. The main products in 
this region are ^SSi, ^28, and ^''Ni. CO138E50 pro- 
duces 0.2Mo of ^^Ni in this region. Other important 
species such as ^^Fe, ^^Co, and ^^Mn (decaying into 
^^Cr, ^^Mn, and ^^V, respectively) are synthesized 
more abundantly in the more energetic explosions. 

(4) For the larger explosion energy, oxygen burning takes 
place in more extended, lower density regions. O, C, 
and Al are burned more efficiently in these cases, 
and the abundances of the elements in the ejecta are 
smaller, while a larger amount of ash products such 
as Si, S and Ar is synthesized by oxygen burning. 

6. CONCLUSIONS AND DISCUSSION 

In this paper, we have presented a model for SN 1998bw 
which is in better agreement with the early observations 

^After submitting this paper, Sollerman et al.'s paper (2000) was published. They found that ^''Ni of ^ 0.3 - 0.9Mq was necessary to power 
the late light curve based on the model C0138 (IMN98), which is consistent with our results. 



than the previous model in IMN98. Models with differ- 
ent i?K yield different synthetic spectra, and by compar- 
ing with the observed early-time spectra of SN 1998bw and 
trying to fit the very broad absorption features, we selected 
model CO138E50 with Mej = IOMq and Sk = 5 x 10^^ gj-g 
as the best match to the early data. The large value of Ek 
quahfics SN 1998bw as 'the' Type Ic Hypernova. (SN Iln 
1997cy may be called a "Type Iln Hypernova" ; Germany 
et al. 2000; Turatto et al. 2000). The mass of the progeni- 
tor C-l-0 star is I3.8M0, corresponding to a main sequence 
mass of ~ 4OM0. All models require M{^^m)r^ OAMq to 
power the bright light curve peak ^. This is about an or- 
der of magnitude larger than in typical core-collapse SNe. 
The compact remnant is probably a black hole, because 
its mass exceeds ^ 3Mq as constrained from the mass of 
56Ni. 

Although the early fight curve of SN1998bw {t < 50 
days) is reproduced well by the most energetic model 
CO138E50, the observed tail declines more slowly than 
this model does {^). The lower energy model CO138E30 
is in better agreement with observations at t < 100 days 
but it still declines is too fast. Model C0138E7 has a 
slower tail and can reproduce the observed light curve tail 
if a smaller ^^Ni mass of 0.28 M© is adopted. This suggests 
that there might be a high density core with low velocities 
in SN1998bw where the 7-rays deposit efficiently. 

Such a dense core could be formed in spherically sym- 
metric models if the exploding star had a massive He en- 
velope so that a strong reverse shock was formed at the 
C+O/He interface and largely decelerated the inner core 
as found in SN lb models (Hachisu et al. 1991, 1994). 
However, the apparent absence of the He I lines in the 
optical spectra in SN 1998bw (Patat et al. 2000) is not 
consistent with the presence of such a massive He enve- 
lope. Also late time energy input from ^^Co bubbles might 
create a non-homologous structure. However, the energy 
from radioactive decays is much smaller than the large ex- 
plosion energy of SN 1998bw {^ 10^^ ergs), so that its 
effect is negligible. 

We suggest that the peculiar density distribution is the 
result of a non-spherically symmetric explosion. If the out- 
burst in SN1998bw took the form of a prolate spheroid, for 
example, the explosive shock was probably strong along 
the long axis, ejecting material with large velocities and 
producing abundant ^^Ni, which might have caused the 
early bright light curve. In directions away from the long 
axis, on the other hand, oxygen would not be consumed, 
and the density could be high enough for 7-rays to be 
trapped even at advanced phases, thus giving rise to the 
slowly declining tail. These features can be seen in the hy- 
drodynamical models of jet-like explosions (Maeda et al. 
2000; MacFadyen & Woosley 1999; Khokhlov et al. 1999; 
Nagataki et al. 1997). The fact that the late light curve 
is fitted by model C0138E7 with 0.28 Mq of ^^Ni, while 
the early light curve requires as much as 0.4 Mq, could 
also support this suggestion, and point at the geometrical 
ratio between the extension of the high- and low-density 
regions. 

There are several observations which support the above 
non-spherical explosion scenario. The observed polariza- 
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tion (- 0.5%) of the early optical hght (Patat et al. 2000; 
IMN98) suggests that the ejecta of SN 1998bw is aspher- 
ical. The observed abnormal distribution of elements in 
velocity space as seen in the emission line profiles in the 
nebular phase, with significant amounts of O being located 
at lower velocity than Fe (Patat et al. 2000; Danziger et 
al. 1999; Nomoto et al. 2000) can better be explained if 
the explosion is non-spherically symmetric (Maeda et al. 
2000). 

The possible connection between SN 1998bw and 
GRB980425 also supports the conjecture that SN 1998bw 
was aspherical. The energy of the photons produced by 
synchrotron emission at the relativistic shock is approxi- 
mately given by hiy ~ 160 keV (r/lOO)-*/!^^ (Piran 1999), 
where F is the Lorentz factor of the shock and ni (cm^'^) 
is the density of the interstellar matter. In order to pro- 
duce an observable GRB, F should be as large as F '^ 100. 
However, even the most energetic model (CO138E50) has 
only a very small mass of relativistic ejecta (~ lO~^*'A/0 
with F > 100), although relativistic hydrodynamical cal- 
culations are necessary to obtain the accurate mass. Such 
a small amount of relativistic material, which is consistent 
with previous estimates (IMN98; Woosley et al. 1999), is 
not enough to produce GRB980425 in a spherically sym- 



metric model. However, if the explosion is axi-symmetric, 
for instance, the energy can be carried by only a small 
fraction of the material, which might then attain a large 
Lorentz factor. 

Regarding asphericity, we note that Danziger et al. 
(1999) and Mazzali et al. (2000b; see also Nomoto et al. 
2000a) estimated a ^'^Ni mass of 0.35 - 0.65 Mq from the 
nebular lines of Fe of SN 1998bw. This estimate does not 
depend much on the asphericity, and is in good agreement 
with the ^^Ni mass of the spherical models C0138. On the 
other hand, Hoflich et al. (1999) suggested that the ^^Ni 
mass can be as small as 0.2 M© if aspherical effects are 
large. Because the difference between these results is not 
so large, aspherical effects might be modest in SN 1998bw. 
If the ^^Ni mass could be determined more accurately from 
the late observations, it would provide a good measure of 
the degree of the asphericity. 

We would like to thank Ferdinando Patat and Enrico 
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Ministry of Education, Science, Culture, and Sports. 
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Fig. 1. — Density distributions against the velocity of homologously expanding ejecta (above) and the velocity profiles ageinst the enclosed 
mass (below) for CO138E50 (soUd line), CO138E30 (long-dashed line), C0138E7 (short-dashed line), and C0138E1 (dash-dotted line) at t = 
250 sec. 
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Fig. 2.— The light curves of CO138E50 (Ek = 5 x 10^2 erg; sohd hne) and CO138E30 {Ek = 3 x lO^^ ^.g; long-dashed line) compared 
with the bolonietric light curve of SN1998bw (Patat et al. 2000). A distance modulus of /i = 32.89 mag and Ay = 0.05 are adopted. The 
dotted line indicates the energy deposited by positrons for CO138E50. The HST observation at day 778 (Fynbo et al. 2000) is shown by 
assuming negligible bolometric correction. 
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Fig. 3.— The light curves of C0138E7 {Ek = 7 x 10^^ erg; dashed Une) and C0138E1 {Ek = 1 x 10^^ erg: dash-dotted Hne) compared 
with the bolometric Ught curve of SN1998bw (Patat et al. 2000). Also shown is the light curve of modified C0138E7 with smaller ^®Ni mass 
of O.28M0 (dotted line). The HST observation at day 778 (Fynbo ct al. 2000) is shown by assuming negligible bolometric correction. 
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Fig. 4.— Photospheric velocities of CO138E50 (solid line), CO138E30 (long-dashed line), C0138E7 (short-dashed line), and C0138E1 
(dash-dotted line) compared with the observations of SN1998bw. Shown are the photospheric velocities obtained from spectral models (filled 
circles, this paper), the observed velocity of the Si II 634.7, 637.1 nm lines measured in the spectra at the absorption core (open circles; Patat 
et al. 2000), and that of the Ca II H + K doublet measured in the spectrum of May 23 (square; Patat & Piemonte 1998). 
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Fig. 5. — The eraly light curves of model CO138E50 with modified ®^Ni distribution (sec text) compared with the bolometric light curve 
of SN1998bw (Patat et al. 2000). 
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Fig. 6. — Three ESO spectra of SN 1998bw near maximum (bold lines) are compared with synthetic spectra obtained using model CO138E50 
(continuous lines) and with spectra obtained for a modified density distribution as described in the text (dotted lines). 
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Fig. 7. — The isotopic composition of cjecta of hypernovae (Ek = 5x10^^ ergs; top left), (_Ek = 3x10^^ ergs; top right), (_Ek = 7x10®^ ergs; 
bottom left), and usual supernovae {Ej^ = 1 x 10^^ erg; bottom right). Only the dominant species are plotted. The explosive nucleosynthesis 
is calculated using a detailed nuclear reaction network including a total of 211 isotopes up to ^^Ge. 
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Table 1 

Parameters of the C+O star models. 



model 


C 







Ne 




Mg 


Si 


S 


Ca 


Ti 


Fe 


Ni 


C0138E1 


0.10 




8.9 


0.66 




0.51 


0.52 


0.19 


0.023 


0.0007 


0.44 


0.065 


C0138E7 


0.09 




8.5 


0.58 




0.47 


0.52 


0.19 


0.029 


0.0010 


0.44 


0.036 


CO138E30 


0.07 




7.6 


0.44 




0.38 


0.74 


0.35 


0.054 


0.0009 


0.45 


0.023 


CO138E50 


0.06 




7.1 


0.37 




0.34 


0.83 


0.41 


0.066 


0.0007 


0.46 


0.018 


model 


44rpi 




""Ni 


f'Ni 


















C0138E1 


3.4x10- 


-4 


0.40 


1.6x10" 


-2 
















C0138E7 


4.5x10" 


-4 


0.40 


1.4x10" 


-2 
















CO138E30 


2.3x10" 


-4 


0.40 


1.2x10" 


-2 
















CO138E50 


5.5x10" 


-5 


0.40 


1.1x10" 
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Table 2 
Yields of hypernova and supernova models (M©) 



i2c 


6.45E-02 


i3c 


2.57E-08 


i4N 


1.34E-07 


ibN 


4.08E-08 


ibQ 


7.07 


1^0 


1.76E-07 


180 


2.45E-08 


19p 


1.12E-08 


20 Ne 


3.72E-01 


21 Ne 


l.lOE-03 


^^Ne 


5.49E-04 


23Na 


1.18E-02 


^4Mg 


2.57E-01 


25Mg 


2.76E-02 


26Mg 


5.78E-02 


^Ul 


4.99E-02 


28Si 


7.27E-01 


29Si 


4.30E-02 


30Si 


5.70E-02 


31p 


6.93E-03 


32g 


3.88E-01 


33s 


1.48E-03 


34g 


2.46E-02 


36s 


2.12E-05 


35 ci 


5.17E-04 


37C1 


1.39E-04 


36 Ar 


7.10E-02 


38 Ar 


l.OOE-02 


40 Ar 


2.85E-07 


39k 


4.26E-04 


4lK 


3.21E-05 


4"Ca 


6.61E-02 


42Ca 


2.83E-04 


43Ca 


4.60E-07 


44Ca 


5.53E-05 


4«Ca 


5.53E-10 


48Ca 


7.90E-14 


45 So 


1.02E-06 


46Ti 


1.09E-04 


47Ti 


2.96E-06 


4«Ti 


5.87E-04 


49Ti 


4.90E-05 


50Ti 


4.15E-09 


soy 


l.OOE-08 


Sly 


8.39E-05 


50Cr 


5.11E-04 


52 Cr 


9.53E-03 


53Cr 


1.16E-03 


54Cr 


3.07E-07 


55 Mn 


4.32E-03 


54pg 


4.92E-02 


56pg 


4.01E-01 


'"Fc 


1.07E-02 


58pg 


1.04E-07 


59C0 


1.72E-04 


58Ni 


1.36E-02 


60Ni 


3.54E-03 


6lNi 


8.69E-05 


62Ni 


5.91E-04 


64Ni 


3.10E-13 


63CU 


3.47E-07 


s^Cu 


2.36E-07 


64Zn 


6.49E-06 


66Zll 


3.91E-06 


e^Zn 


3.03E-09 


68Zn 


1.93E-09 



















Table 3 
Detaied yields of model CO138E50 (Mq) 



